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Abstract 

We review some aspects of the fermionic interpretation of the two-dimensional Ising 
model. The use is made of the notion of the integral over the anticommuting Grassmann 
variables. For simple and more complicated 2D Ising lattices, the partition function can 
be expressed as a fermionic Gaussian integral. Equivalently, the 2D Ising model can 
be reformulated as a free-fermion theory on a lattice. For regular lattices, the analytic 
solution then readily follows by passing to the momentum space for fermions. We also 
comment on the effective field-theoretical (continuum- limit) fermionic formulations for 
the 2D Ising models near the critical point. Key words: Ising model, anticommuting 
Grassmann variables, critical point. PACS number(s): 64.60. Cn; 75.10.-b; 05.30.Fk.[] 

1 Introduction. 

The two-dimensional (2D) Ising model in a zero magnetic field was first solved by Onsager 
in his celebrated paper [1]. This remarkable solution have already contributed much to our 
understanding of the nature of the second-order phase transitions in magnets [2-6]. The 
original algebraic method used by Onsager, however, was enough complicated. The modern 
approaches to the 2D Ising model are based merely on the fermionic interpretation of this 
model [7-23]. The fermionic structures in the 2D Ising model were first recognized within 
the transfer-matrix and combinatorial considerations [7,8]. It was realized later on that the 
notion of the integral over anticommuting Grassmann variables due to Berezin [9,10] is a 
powerful tool to analyze the 2D Ising models [10-23]. Berezin himself was the first to apply 
the anticommuting variables to evaluate the partition function of the 2D Ising model on 
rectangular lattice [10]. His method was based on combinatorial considerations [10]. The 
fermionic approaches to the 2D Ising model have been developed later on by many authors 
[11-23]. In the forthcoming discussion we will merely follow a simple fermionic interpretation 
of the 2D Ising model introduced in [17-19]. The approach is based on integration over 
the anticommuting Grassmann variables and the mirror-ordered factorization principle for 

1 Published in: J. Phys. Studies (Ukr), Vol. 1, No. 4 (1997) p. 554-563. 
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the 2DIM density matrix. The method do not involve the traditional transfer-matrix or 
combinatorial considerations. Schematically, we have: 

Q= SpQ(<r)-> Sp Q{a\a) -> Sp Q (a) = Q. (1) 

(a) (a | a) (a) 

Here we start with the original spin partition function, Q, and introduce, in a special way, 
a set of new anticommuting Grassmann variables (a) thus passing to a mixed (a \ a) rep- 
resentation for Q. Eliminating the Ising spin variables in this mixed (a \ a) representation, 
we obtain a purely fermionic expression for the same partition function Q. The partition 
function appears, in final form, as a Gaussian fermionic integral. Equivalently, the 2D Ising 
model is represented as a free-fermion field theory on a lattice [17-19]. For the homogeneous 
(translationally invariant) lattices the analytic solution for the partition function and free 
energy then readily follows by passing to the momentum space for fermions by means of 
Fourier substitution. In particular, this gives a few line derivation of the Onsager result 
[17]. The exact solution for a finite lattice on a torus (periodic boundary conditions in 
both directions) also follows within suitable modification of the fermionization procedure 
[18]. The decomposition of the partition function on a torus into a characteristic sum of 
four fermionic integrals follows here from a general rule of transposition of two Grassmann 
variable functions derived in Appendix A of [18]. The Ising models settled on the lattices 
with complicated local structures have being analyzed by factorization method within the 
spin-polynomial interpretation of the problem [19]. We will pay more attention to this last 
case in sections 3-7. For application of the results of fermionic analysis [19] in studies of 
regularly diluted Ising ferromagnets also see [25,26]. It is also important that the method 
works in fact for the most general inhomogeneous distribution of the coupling parameters 
over the lattice bonds [17]. This may be of interest in studies of the disordered Ising models 
[20-23]. The Gaussian fermionic representations has been recently constructed as well for 
the inhomogeneous two-dimensional dimer problems [24]. 

In the straightforward variant of the method [17] we introduce Grassmann variables al- 
ready at the first stages in order to decouple the local bond Boltzmann weights into separable 
factors called Grassmann factors (GFs). We then combine GFs with the same spin variables 
into separable groups all over the lattice and then sum over spin states in each group indepen- 
dently, thus passing to a purely fermionic representation for Q. In general, the factorization 
idea we apply resembles the idea of insertion of the Dirac unity, £ | a) (a | = 1, in quantum 
mechanics. However, the situation in our case is more complicated since we deal here with 
Grassmann factors which are neither commuting nor anticommuting, in general, with each 
other. The key point of all the construction is then the mirror-ordering procedure for the 
global products of GFs [17]. This ordering procedure enables us to keep nearby few relevant 
factors with the same spin variable in the process of fermionization, which is the necessary 
condition in order we can actually eliminate spin degrees of freedom. The ordering procedure 
for global product of GFs is the most important part of the solution [17-19]. It might be 
worthwhile also mentioning that the lacking of a suitable ordering procedure for GFs in the 
3D case is a serious obstacle to deal with the 3D Ising model within fermionic factorization 
method. The same holds true for the 2D Ising model in a non-zero magnetic field. 

If we start with factorization of the local bond Boltzmann weights, as we do for the stan- 
dard rectangular lattice in [17,18], we obtain the fermionic representation for Q with four 
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variables per site. Starting by factorization of the bond weights, the number of fermionic 
variables per site will be even larger for more complicated lattices. An important modifi- 
cation of the method was introduced in [19], where we start with factorization of the local 
cell weights presented by three-spin polynomials. The spin-polynomial interpretation for 
the 2D Ising models arises if we multiply few local weights forming elementary cell in Q. 
Effectively, this reduces the number of fluctuating degrees of freedom already on the level 
of the spin-variable formulation of the problem, prior to fermionization. Respectively, the 
Gaussian integral for Q appears here with only two fermionic variables per site, see ([15]) be- 
low, which provides essential simplifications in the forthcoming analysis [19]. For a set of 2D 
Ising models, including rectangular, triangular and hexagonal Ising lattices as the simplest 
cases, this gives analytic solution in terms of the parameters of the three-spin polynomial 
characterizing elementary cell. The exact lattice integral ( 15]) is also a suitable starting point 
to pass to the effective continuum-limit formulation near T c . The continuum-limit formu- 
lation corresponds to the low-momentum sector of the exact lattice theory responsible for 
the critical-point singularities and the large-distance behaviour of correlations near T c . The 
resulting theory is the massive 2D Majorana theory, with mass vanishing at T c . By doubling 
of fermions in Majorana representation, we can pass as well to the 2D Dirac field theory of 
charged fermions. Below we comment on some aspects of the fermionic interpretation of the 
2D Ising model in more detail. The fragments of present discussion have been reported also 
in [27,28]. In the next section we outline shortly the basic rules of fermionic integration. 



2 Grassmann variables. 

We remember that Grassmann variables are the purely anticommuting fermionic symbols. 
Given a set of Grassmann variables a±, ci2, a/v, we have a^aj + ajcii = 0, a| = 0. The 
Beresin's rules of integration for one variable are [9]: 

J ddj ■ cij — 1 , J dctj -1 = 0. (2) 

In the multidimensional integral, the differentials dai, da-i, ... , da^ are again anticommuting 
with each other and with the variables. The basic relations of the grassmannian analysis 
concern the Gaussian fermionic integrals [9,10]. The Gaussian integral of the first kind is 
related to the determinant: 

N I N N \ 

/ Y\ da*da,j exp ciiA^a* = det A , (3) 

3=1 V=lj=l / 

where {aj,a*} is a set of completely anticommuting Grassmann variables, the matrix in 
the exponential is arbitrary. The fermionic exponential here is assumed in the sense of 
its series expansion. The series terminates at some stage due to the property af = 0. The 
correspondent finite polynomial for the exponential in (|D can be obtained also by multiplying 
the elementary factors exp (ajA^a*) = 1 + ^A^a*. The appearance of determinant in (3) is 
due to the known interrelations between fermionic algebra and determinant combinatorics. 
By convention, the variables and a* can be considered as complex conjugated fermionic 
fields, otherwise these are simply independent variables. The Gaussian integral of the second 
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kind, for real fermionic fields, is related to the Pfaffian of the associated skew-symmetric 
matrix: 

r (1 N N \ 

/ dajsi ... da2da\ exp - EE a *^k? a i = Pfeff^j = ■ (4) 

V 2 »=ij=i / 

The Pfaffian is some combinatorial polynomial in elements known in mathematics for a 
long time. Under another name, the Pfaffian is also well known in physics since the Pfaffian 
combinatorics is identical with that of the fermionic version of Wick's theorem. For any 
skew-symmetric matrix (Aij = —Aji) we have: 

deti= (Pfaffi) 2 . (5) 

This identity can be most easily proved just in terms of the fermionic integrals like (H) 
and (^). The fermionic averages (Green's functions) correspondent to the integrals (^) and 
(fD can be defined in a natural way. The linear superpositions of Grassmann variables are 
again Grassmann variables and it is possible to make a linear change of variables in the 
fermionic integrals. As compared with the rules of commuting analysis, the only difference 
is that Jacobian will now appear in the inverse power [9,10]. [It might be also interesting to 
compare (|3|) with the analogous Gaussian integral over the commuting complex- variable fields 
in common analysis]. In physical applications, the new fermionic variables of integration are 
often introduced by Fourier substitution (transformation to the momentum space). 



3 The spin-polynomial interpretation. 

In order to construct the fermionic representation for Q with the minimal number of fermionic 
variables per site, we start with the spin-polynomial formulation of the problem [19]. This 
can be most simply illustrated by an example of the Ising model on a triangular lattice. The 
triangular lattice can be viewed as a rectangular one with an additional diagonal interaction 
introduced in each rectangular cell. Let lattice sites on the corresponding rectangular net 
be marked by integer Cartesian coordinates mn , with m — 1,2, L, n = 1, 2, L running 
in horizontal and vertical directions, respectively. Here L is the lattice length, N = L 2 is 
the total number of sites and spins in a lattice. At final stages we assume L 2 —>■ oo. With 
each lattice site we associate the Ising spin variable a mn = ±1. The hamiltonian of the Ising 
model on a triangular lattice is then given as follows: 

— (3 H (a) = ^ [ b\ (Tmn&m+ln + &2 ffm+ln&m+ln+l + ^3 &mn a 'm+ln+l } > P = 1/ kT , (6) 

mn 

where b a = J a /kT are the dimensionless coupling constants, J a are the magnetic exchange 
energies, and kT is the temperature in energy units. To be definite, in what follows we will 
merely keep in mind the ferromagnetic case, b a > 0, though this restriction is not essential 
until we pass to the continuum limit formulations. The partition function and the free energy 
per site are: 

z = E«" w - -0fz= i im 4 lnZ ' ( 7 ) 

O) 
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where the sum is taken over 2 N spin configurations provided by a mn = ±1 at each site. 

Noting the identity for the typical bond weight: e bcrcr = cosh b + sinh b ■ aa', which 
readily follows from (erer') 2 = +1, we find: 

Z = (2 cosh bi cosh b 2 cosh b 3 ) N Q , N = L 2 -> 00 , (8) 

where Q is the reduced partition function: 

q= s P { n ( 1 + *i )(l+*2 C"m+ In ""m+ ln+ 1 O mn < ~> m+ln+1 )}, (9) 

( CT ) mn 

where i Q = tanh 6 Q , 6 a = J a /kT, and S^o-) is the properly normalized spin averaging [19]: 

Sp (...) = n s p (-) = 2^E(-). sp (...) = i £ (•••), (10) 

( CT ) mn. (ffmn) (omn) ^ CTmn=±1 

the local averaging is normalized here in such a way that Sp (1) = 1, also note that Sp (cx mn ) = 
0. The reduced partition function Q will be the main subject of our interest. 

In order to pass to the spin-polynomial interpretation, we have to multiply the three 
bond weights forming a triangular cell in Q. Let us introduce the local enumeration of sites 
at our rectangular net: 

( C"l I C"2 I CT3 I C4 ) *-* ( fmn | a m+ln | 0" m+ln+1 | O mn+1 ) • ( 1 1 ) 

The elementary cell in (§) and (^|) is formed by a triangle of spins (eri, cr 2 , cr 3 ) mn . Noting 
that <jj = ±1, and hence a 2 = 1, we can multiply the weights in @ making use of the 
properties like o\(J2 • o"20"3 = o"iO"3, etc. For the product of the three bond weights from (f|) 
we then obtain the three-spin polynomial as is given in fll2"D below. This polynomial can 
be considered as the effective Boltzmann weight of a triangular elementary cell taken as a 
whole: 

-Pl23(cr) = (1 + h 0i<7 2 ) (1 + h 0-20-3) (1 + *3 0-10-3) 

= (1 + ^^3) + (*! + t 2 t 3 ) oia-2 + (t 2 + tit 3 ) 0203 + (t 3 + *i* 2 ) 0103 (12) 
= a + «i 0!cr 2 + « 2 o 2 cr 3 + a 3 a t a 3 . 

In the last line we write the three-spin polynomial in general notation, for the given case of 
triangular lattice we have: 

a = l + t 1 t 2 t 3 , a 1 = ti + t 2 t 3 , a 2 = t 2 + t 1 t 3 , a 3 = t 3 + t 1 t 2 , (13) 

while the case of rectangular lattice follows with t 3 = 0, which corresponds to b 3 = in the 
hamiltonian For rectangular lattice we have: = 1, a\ = ti, a 2 = t 2 , a 3 = t\t 2 . The 
hexagonal and more complicated lattices can be characterized by the three-spin polynomials 
like (|12|) as well, with their own sets of the a-parameters. Thus we come to the three-spin 
polynomial partition function [19]: 

Q= Sp [J (a + ot\ <J mn a m+ i n + ol 2 a m+ i n a m+ i n+ i + a 3 a mn a m+ i n+ i 

) • (14) 

(°-) mn 

In what follows we assume ao, ai,a 2 , a 3 to be arbitrary parameters. Particular Ising lattices 
are specified by the choice of the a -parameters [19,25,26]. 
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4 Free-fermion representation and analytic 
results (lattice case). 

The three-spin polynomial partition function fll4l) can be transformed into a fermionic Gaus- 
sian integral following ([I]). The starting point is factorization of the local polynomial weights 
in fll4[) . We need only two variables per site in the process of fermionization. The partition 
function then appears as the following Gaussian fermionic integral [19| : 

Q = Yl d c mndc mn exp ^ I ( a c mnC* mn ~ «1 Cmn C m-ln ~ a 2 CmnC-1 
mn mn 

0^3 Cmn.C m _i n _i) OL\ C mn C m —i n Q.2 C-mn^mn—X ] , (15) 

where c mn , c* mn are the totally anticommuting Grassmann variables, two per site. The free- 
fermion representation ([To]) is exact. The integral fll5|) is equivalent to the original partition 
function ([TJ]) up to the boundary effects, which are inessential as L 2 — > oo. The explicit 
evaluation of the integral fll5|) can be performed by passing to the momentum space by 
means of Fourier substitution: 

i i 

r = —X^ r P i2 r( m P+ nc l) r * = — r * p ~ i 2 r( m P+ n l) flfil 
u mn j / j u »q J i ^mn r VI ' \ ±KJ J 

Here c , c* q are the new fermionic variables of integration. The integral ( |T5D now appears 
in the form: 



Q = [[dc* q dc pq exp 2^ [ c pq c* q [a - a x e ! ~ - a 2 e l ~ - a 3 e* ~ 



L (p+g) 

/»/ /"/ 

- ai c w c L _ pL _ g e 1 -r - a 2 c* L _ pL _ q c* q e 1 "z 



(17) 



Note that the fermionic measure transforms in a trivial way by passing from (|15| ) to ( pTD . This 
is because the Jacobian of the substitution (|16D is unity, as it follows from the orthogonality 
of the Fourier eigenfunctions. 

It is easy to see that the integral fllTD decouples into a product of simple low-dimensional 
integrals. Since the variables with momenta pq and L — pL — q interact in ([T7|), we have to 
single out in the fermionic action in ( PD the combined (pq \ L — pL — q) term. Let S pq be 
the term already given explicitly in the pq sum above. It then follows that the integral (|17|) 
is the product of the following independent factors: 

Q\ q = J D vg ex P ( S vg + S L - P L-g) , D pq = dc* L _ pL _ q dc L _ pL _ q dc* pq dc pq . (18) 

The elementary integral ([IB]) can be readily evaluated by elementary tools (for instance, 
simply using the definitions (Q)) and its value is given as a factor in the product ( fH]| ) below. 
By comparing the fermionic measures in (|T7|) and (|18D, it follows that the partition function 
itself, Q, arises if we multiply the factors fllq) over only one half of the points in the mo- 
mentum space. That is, we have to multiply the factors ([18]) in such a way that if the given 
mode pq is already included into the product, then the conjugated mode L—pL — q is not to 
be included, and vice versa. Respectively, the total product of factors flTB"|) over the complete 
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set of momentum modes, < p, q < L — 1, will yield the squared partition function. Thus, 
for the partition function fll5| ) we find [19]: 

L-1L-1 

Q 2 = II II [ ( a + a l + a 2 + a D - 2 ( a a l - a 2 a 3) COS —j- 

p=0 q=0 L 

2nq 2n(p + q) 
— 2 (a a 2 — ce^z) cos — 2 {a a 3 — a 1 a 2 ) cos J . (19) 

Lj Lj 

Respectively, the free energy per site is: 

-Pf Q = Km (^lnQ) 

1 f J dp dq 2 2 2 , 2^ 



2 J J zn Zn ' 



— 2 (a a;2 — a i a ^) cosg — 2 (a a 3 — a 1 a 2 ) cos(p + g) ] . (20) 

In particular, the exact solutions for the rectangular, triangular, and hexagonal Ising lattices 
follow from ( |2"0"| ) under correspondent specifications of the parameters a. [19]. 

Here we consider for illustration the simplest case of the standard rectangular lattice. 
For this lattice, we have a , at, a 2 , a 3 <-> 1, ti,t 2 , t\t 2 . From (120), the free energy is: 



2?r 2tt 



Mq I r ect = ^ / / ^ [ (1 + t\) (1 + # - 2*! (1 - t\) cosp 







— 2t 2 (1 — tx) cos q ] . (21) 



The free energy (^TJ) is associated with the reduced partition function, Q, while the true 
partition function is Z = (2c\C2) N Q, see (|j) and @ with £ 3 = 0. Also we remember the 
identities like c 2 a {l + t 2 a ) = cosh26 Q , 2t a c 2 a = sinh26 Q , etc, with b a = J a /kT. From (pl|), the 
true free energy per site then appears in the form: 

1 f f dp dq 



(3fz | rect = ln2+ - / / In [ cosh26i cosh26 2 — sinh26i cosp 

2 J J 2ix 2ix 



o o 



sinh26 2 cosg] . (22) 



This is the famous Onsager solution for the free energy of the 2D Ising model on a rectangular 
lattice [1] . An interesting comment on the structure of this solution follow immediately after 
equation (108) in [1]. 

As it follows from (|21|), for rectangular lattice the critical point is fixed by condition: 



1 - t x - t 2 - t x t 2 = , t a = tanh( J a /kT) , (23) 
equivalently, this condition can be written in the form: 
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sinh(2Ji/fcT) sinh(2J 2 /fcT) 



(24) 



which rather corresponds to the solution in the form of (22). In both formulations the 
ferromagnetic case is assumed, 61,62 > 0, equivalently, ti,t 2 > 0. The free energy is known, 
the specific heat can be obtained by differentiation with respect to the temperature: C/k = 
f3 2 (<9 2 (— flf)/d(3 2 ), where C/k is the dimensionless specific heat, (3 = 1/kT, and k is the 
Boltzmann constant. The singularity in the specific heat appears to be logarithmic near T c : 
C/k ~ I In |t||, ' 



(T — T c )/T c —> 0, as it can be deduced ^from fl2ll) and/or (|22l). Also see 



the discussion in the next section. 



5 The symmetries and the critical point. 



The symmetries provided by the above exact solution ([T9|) and (|20|) , and the closely related 
question on the location of the critical point, have an interesting interpretation within the 
spin-polynomial language [19]. Making use of the local enumeration of sites in the elemen- 
tary cell, see (111]); the cell weight in the density matrix in (|14]) is given as the three-spin 
polynomial: 

^123 (0-) = a + ai ai a 2 + a 2 a 2 a 3 + a 3 a x a z . (25) 
It appears also to be useful to introduce the associated three-spin polynomial: 

-F123 (0) = a - a x o x o 2 - a 2 a 2 a 3 - a 3 ai<j 3 , (26) 
and we note the following interesting identity [19]: 

(-^123(0)) 2 = ("o - Oil o x o 2 - a 2 o- 2 03 - a 3 a x a 3 ) 2 
= (q!q + a\ + a\ + ctg) — 2 (a ai — a 2 a 3 ) a 1 a 2 
- 2 (a a 2 - a x a 3 ) a 2 a 3 - 2 (a a 3 - a x a 2 ) a x a 3 . (27) 

It is seen that the combinations of the a-parameters occurring in (-F123) 2 are exactly the 
same as in the momentum modes Q 2 (p \ q) in (p0|): 

Q 2 (p\q) = (chq + a\ + a\ + a\) — 2 (a a 1 — a 2 a 3 ) cosp 
—2 (a a 2 — a 1 Q!3) cos q — 2 (a Q a 3 — a x a 2 ) cos(p + q) , (28) 

where < p, q < 2ir, the limit L 2 — ► 00 is already assumed. 

It appears that the following combinations of the a-parameters play important role in 
discussing the symmetries and the critical point [19]: 



a o — \ ( a o + oil + a 2 + a 3 ) , a = \ (a — a x — a 2 — a 3 ) , 

Oi\ — \ (a + a x — a 2 — a 3 ) , a 1 = \ (a — «i + a 2 + a 3 ) , 

a* 2 = \ (a — «i + a 2 — a 3 ) , a 2 = \ (a + «i — a 2 + a 3 ) , 

«3 = 5 («o — Oil — ol 2 + a 3 ) , a 3 = \ (a + Oil + a 2 — a 3 ) . 
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The parameters a* and a are in fact the eigenvalues of the polynomials \P122, and 5^123, 
respectively. By the "eigenvalues" we mean the four numbers which takes the polynomial as 
the spin variables run over their permissible values ±1. We note also the following important 
identity [19]: 

OiQOna2a3 = a^ala^al — aoaic^cts . (30) 



There are some evident symmetries in the solution fl2~0"|). For instance, the free energy (20) is 
a symmetric function with respect to arbitrary permutations of ao, cx\, «2, a 3 parameters. We 
can as well change the signs of any two of them, with — /3/q unchanged. There is also a less 
evident hidden symmetry in the solution, corresponding to the Kramers- Wannier duality in 
the case of the standard rectangular lattice. Namely, the partition function Q{a} is invariant 
under the transformation a, <-> a*. This symmetry in fact holds already for the parameters 



of the separable fermionic modes (28) , and can be proved making use of (|27|), for further 
details see [19]. 

In order to establish the possible critical points, we have to look for zeroes of Q 2 (p \ q) 
momentum modes. As it can be guessed already from the analogy between ( p8|) and (|27|) , 
there are four exceptional modes with (p \ q) — (0 | 0), (0 1 7r), (-7T | 0), (ir \ ir). For these modes 
we have, respectively, Q 2 {p\q) = (2a ) 2 , (2<5i) 2 , (2a 2 ) 2 , (2a 3 ) 2 . Remember that the 
parameters and hence aj are some functions of temperature. Thus, if at some temperature 
one of the above momentum modes vanishes, we fall at the point of phase transition. It can 
be shown that all other modes Q 2 {p\q) are definitely positive at all temperatures, there 
are no other critical points for physical values of the a-parameters. The possible criticality 
conditions can be combined into one equation: 

OioCii 0(2 &3 = 0. (31) 

It can be shown also that if all bond interactions are ferromagnetic then a%, a.2, a.3 never 
become zero, and the critical point can only be associated with «o = 0. Thus, for the purely 
ferromagnetic interactions the unique critical pont is definite by the equation: 

2 «o = cto ~ a i ~ a 2 — a 3 — . (32) 

The criticality conditions with a.^ — 0, k — 1, 2, 3, can only be realized when the antiferro- 
magnetic interactions are involved in the hamiltonian. 

If the critical point is associated with 6tj = 0, then the singular part of the free energy 
(p0|) near T c is given as follows [19]: 



(2a j) 2 , const 

= In — — - 

16irJ(a aia2as) c V' a ii 



PfQ I singular _ \ „ f, " ^ ^ n /o^. \2 ' (^) 



which implies that, near T c , in the order of magnitude — /3Jq ~ r 2 lnr 2 , where r ~ \T — T c \. 
The specific heat have thus the log-type singularity, C ~ | lnr | as T — > T c (Onsager, 1944). 

It is seen that the eigenvalues ( p5|) play important role, but it is not so clear, in fact, 
how the role of the polynomial F123 can be understood in a less formal way, at the level of 
the original spin-variable formulation of the problem, prior to the analytic solution. The 
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special role which the parameters (|29| ) play in the inherent structure of the 2D Ising model 
is confirmed by the expression for the spontaneous magnetization. The 8-th power of the 



spontaneous magnetization M = (cr mn ) for model (|T1| ) is given by the following very simple 
expression [19]: 

aoaia 2 ct3 aoaia 2 as 

This expression for M 8 holds true when the right hand side varies between and 1, and 
M 8 = otherwise. It is easy to check that the known expressions for the spontaneous 
magnetizations of the rectangular, triangular, and hexagonal lattices follow easily from (j34[) 
as particular cases. From (^) we find M ~ r s as r ~ | T — T c | — >• 0, with the universal 
value of the critical index (3 = 1/8 for the magnetization at the critical isobar. What are the 
hidden reasons for such simple expression for M 8 , this is yet unknown. 



6 Majorana fields. 

The fermionic integral for Q given in (|15|) also appears to be a suitable starting point to 
formulate the continuum limit field theories for the correspondent models near T c p7| . For 



other approaches to the contnuum-limit formulation of the 2D Ising model also see [IR 



22|. Assuming the purely ferromagnetic case, with a ,«i,«2,«3 all positive, we write once 



again the exact lattice action from (|T^) as follows: 



Sic) 



E 



(a - a x - 
+ a 2 

+ «l C, 



a 2 
(c* 



r 

mn mn 



c 



mn—1, 



C 



m—ln) 



+ Ot^Crr, 



(35) 



with cjL 



,* 2 



m -^mn 



•Em— In j 9 n X mn 



X r 



. Let us define the lattice derivatives as follows: d m x 
n ^mn-i ? also note that x mn x m — in— i d m x mn ~t~ n x 



X r 



Introducing the new notation for the fermionic fields, c mn , c* mn 
action (BSD in the form: 



we find the 



$ C0) ~~ E [ TLktymn^mn + ^1 ^Pmn dm ^Pmn + ^2 VWi 9 n 1pm: 
mn _ _ _ 

^Pmn 9 m n Iprnn + "1 i>mn 9 m ^ inn mn I 



(36) 



Ai = QJi + 0:3, A2 = a 2 + «3 j m — ~ a l ~ a 2 — ^3 = 2 «o • 

It is easy to recognize in this still exact lattice action the typical relativistic field-theoretical 
like structure with mass term and kinetic part. Since we assume the ferromagnetic case, the 
criticality condition is 2d = , see (|3~2j). The parameter m = ceo — «i ~ ol 2 — a 3 = 2q 
plays the role of mass in the field-theoretical interpretation. The mass vanishes at the critical 
point, m ~ | T — T c | — > as T — > T c . The ordered phase corresponds to m < . 

Taking in (|36|) the formal limit to the continuum euclidean space, with (mn) — > x = 
(xi\x 2 ), ip mn -> ip(x) = i/j(x 1 \x 2 ), and d m di = d/dx x , d n -> d 2 = d/dx 2 , and 
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neglecting the second-order kinetic term with d\d%, we obtain the Majorana like fermionic 
action of the correspondent continuum theory: 



S 



d 2 x [mip(x)ip(x) + Ai if>(x) d\ if>(x) + X 2 ip(x) d 2 if>(x) 
+ «i i/j(x) di i>(x) + a 2 i>{x) d 2 4>(x) } , 



(37) 



where the parameters Ai, A2, 221 are defined in ([31]). A nonstandard feature in this action are 
the " nondiagonal" kinetic terms ipdiip, ipd 2 %f). These terms can be eliminated by a suitable 
linear transformation of the fields: if) = 7^1 + T]if> 2 , — 701 + ^02 ; where ipi,ip2 are new 
Majorana components. The condition that the nondiagonal like terms will not appear in fl38|) 
is 77 = 7777. Omitting further details, after a suitable transformation of the fields and the axis 
in the d 2 x space, from ( |57D we obtain Majorana action in the canonical form (cf. 



S 



d 2 x [ m *0i (x) i/j 2 (x) + ipi (x) d ip\ (x) - ip 2 (x) d ip 2 (x) ] , 
d = ]- (di + id 2 ) , d = - (d\ — id 2 ) , 2D Majorana , 



(3* 



with the rescaled mass: 



_ «o 
m — - 



a± — a 2 — «3 



1/2 



(2 J (atoatia 2 a 3 ) 



(39) 



The 2D Ising model is presented now as a field theory of free massive two-component Ma- 
jorana fermions in euclidean d 2 x space. The new Majorana fields in this representation, 
■01 , if> 2 1 are the linearly transformed fields if) , if) from (j3~7D . The axis of the d 2 x space are 
also rescaled and rotated as we pass ^from fl37D to (|38f) . Respectively, the mass is rescaled 

a x - a 2 - a 3 ) c = 0. 



according (|39|), here ( ) c means the criticality condition (a 
In matrix notation, the Majorana action (|3l|) becomes: 



S major ~ I d X 



m 



1 
-1 



di+id 2 




Introducing the standard Pauli matrices: 



1 

1 



C 2 





% 



-% 




0-3 





-di + i d 2 



1 







-1 



01 

V'2 



(40) 



(41) 



the action (f40[) can be written also in the form: 
S, 



major 



- J d 2 x^> {x)[m{ia 2 ) +dx{a 3 ) + id 2 (l)} * (s) , 



or in the form: 



x 



01 
^2 



(42) 



5 maj0 r = ^ / d 2 x^{x)[m+d 1 (a 1 ) + d 2 (a 2 )]^{x) , (x) = * (x) (i a 2 ) . (43) 
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Introducing the 2D gamma-matrices in a natural way: 7i = 72 = cr 2 , the action (]43|) 
becomes: 

5 major = ~ / d 2 x^(x)[m +9] , <9 = 7i<9i + 7 2<9 2 , <9 2 = <9 2 + <9 2 , (44) 

which is the 2D Majorana action written in the standard form assumed in relativistic field 
theory. Notice that the conjugated Majorana spinors \P and \& here are not the truly inde- 
pendent fields since they both are built from the same component fields ipi, ip2- The truly 
independent left and right spinors \I/ and ^ appear if pass to the Dirac theory. This is 
commented in next section. 



7 Dirac fields. 

We can pass to the Dirac field theory of charged fermions by doubling the number of fermions 
in the Majorana representation (|38|)-(|44|). On formal level, this corresponds to passing ^from 
Pfaffian like Gaussian integral to the determinantal Gaussian integral according the identity 
|Pfaff A I 2 = det A, see ([5]). To this end, we take the two identical copies S' and 5"' of the 
Majorana action ( [44]) and consider the combined action S'dirac = S' + S". In this combined 
action we introduce the new Dirac fermionic fields \I/ = (ipi \ ?/> 2 ) an d ^ = (V'i* | ^2) by 
substitution: 

¥ = y = -L(&-iV), (45) 

where \1/ ', ^ ", v , ^ are the original Majorana spinors. In terms of the new fields the 
combined action S'dirac — (S + S ) majorana becomes: 

S d irac= Jd 2 x*(x)[m+d]*(x), d = 71 d 1 + 7 2 «9 2 , (46) 



where the gamma matrices are the same as in fl4"4]). If written in components, the substitution 
inverse to ([|5]) takes the form: 

W = -j= (fpi + , » ^1 = -j= (^i - ^2 ) . 
^ = -^(^-■01), iV»2 =-^(^ + ^1)- (47) 



The four components of the Dirac spinors ^ are all independent variables. By convention, 
we can assume ■0*,'0 2 * to be complex conjugated to ipi,ip2- I n this case \l/ is the hermitian 
conjugated to The advantage of the transformation from Majorana to Dirac formulation 
is that we gain in later case the gauge symmetry which can be useful in some cases |[22||. The 



continuum-limit formulations like ( 44[) and (46) captures all relevant features of the exact 



lattice theory at low momenta, or at large space scales, which is only important for the the 
critical-point singularities in thermodynamic functions near T c . For the simplest case of the 
standard rectangular lattice, the continuum limit was used to study the behaviour of the 



2DIM spin correlation functions at large separation in [15, 16]. Another important field of 
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application of the continuum- limit formulation are the Ising models with quenched disorder 
[20-23]. 

In conclusion we note that the singular part of the free energy and specific heat can be 
extracted, in the exact form, from fl44D and/or (fR|). For instance, starting from (f46|) , we 
find: 

-/3/ S in g = \ ^ Q 2 = \ Jd 2 x In Q\x) = ~ J ^ In Q 2 (p) 

= — [ d 2 p\n(m 2 + p 2 ) = -^-m 2 ln + (...) . (48) 

This is equivalent to the equation ( |33"D obtained directly from the exact expression for the 
lattice free energy (^0|), taking into account the definition of the effective mass given in (|39|) . 

8 Conclusions. 

We have discussed some aspects of a simple fermionic interpretation of the 2D Ising models 
in terms of the anticommuting integrals. For any planar 2D Ising model, the partition 
function can be expressed as a fermionic Gaussian integral. The analytic solutions for regular 
lattices then easily follow by transformation to the momentum space. The continuum limit 
field-theoretical formulations for the 2D Ising models also can be deduced readily from the 
exact lattice fermionic integral for the partition function. The differences between particular 
lattices are merely adsorbed, in the field-theoretical limit, in the definition of the effective 
mass. Finally, the free-fermion interpretation for the 2D Ising model shows that this model, 
which in its original formulation is rather discrete combinatorial problem, can be placed, in 
fact, into a common range of the typical models of quantum statistics and solid state physics. 
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